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ABSTRACT 

Finite  element,  elastic  and  plastic  analyses  of  various 
aluminum  panels,  containing  holes  and  notches,  were  conducted 
for  comparison  with  photoelastic  experimental  results.   A 
FORTRAN  IV  program,  ADINA  (Automatic  Dynamic  Incremental 
Nonlinear  Analysis) ,  was  used  for  both  linear  and  nonlinear 
analyses.   Mesh  refinements  were  used  for  each  panel  and  the 
monotonically  convergent  results  were  extrapolated  using 
Richardson's  method.   Stresses  were  locally  smoothed  from 
the  Gauss  integration  points  to  the  nodal  points.   Eight 
noded,  isoparametric  elements  were  used  throughout.   Modi- 
fication to  an  ADINA  preprocessor  program,  also  coded  in 
FORTRAN  IV,  was  made  for  use  with  a  VERSATEC  plotter. 

Comparisons  were  made  to  the  elastic,  analytic  series 
solution  by  Howland  for  a  circular  hole  in  a  finite  strip. 
The  finite  element  results  varied  by  less  than  one  percent 
from  Howland 's  solution.   Handbook  values  for  the  elastic 
stress  concentration  factors  of  the  geometries  investigated 
differ  from  finite  element  results  by  less  than  one  percent 
in  all  cases.   The  photoelastic  works  of  Frocht  were  also 
compared  where  applicable.   Stresses  in  the  plastic  range 
obtained  from  slip-line  theory  for  a  rigid-perfectly-plastic 
material  show  excellent  correlation  to  a  finite  element 
analysis  of  such  a  material.   Comparisons  to  elastic  and 
plastic  experimental  data  were  made  for  the  panels 
analyzed  and  show  good  correlation  to  finite  element  results. 
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I.   INTRODUCTION 

Development  of  on-board  fatigue  monitoring  systems 
for  Naval  aircraft  have  made  it  possible  to  record  exten- 
sive structural  loading  data  in  flight.   The  strain  gages 
used  in  such  a  system  must  be  located  away  from  stress 
concentration  areas  to  prevent  their  fatigue;  however, 
these  areas  are  of  the  greatest  interest  in  analyzing 
and  predicting  fatigue  life  of  the  structure.   Understand- 
ing the  relationship  between  nominal,  far-field  stresses 
and  local  stresses  in  critical  areas  thus  becomes  vitally 
important.   Recent  experimental  investigations  into  the 
effect  of  uniform,  far-field  loads  on  stress  concentration 
areas  have  been  made  at  the  Naval  Postgraduate  School 
(NPS)  using  photoelastic  techniques  [Refs.  1,  2  and  3]. 
These  experiments  involved  loading  7075-T6  aluminum  into 
both  the  elastic  and  plastic  regions,  as  well  as  measure- 
ments of  residual  stresses  resulting  from  plastic  yielding, 

Finite  element  analyses  (FEA)  of  the  aluminum  panels 
used  in  the  experiments  of  Stenstrom  [Ref.  1]  were  con- 
ducted.  The  panels  used  in  the  experiments  of  Engle 
[Ref.  2]  and  Stuart  [Ref.  3]  have  similar  geometry.   The 
finite  element  programs  available  at  NPS  were  surveyed 
and  ADINA  [Ref.  4]  was  chosen  for  its  proven  ability  to 
produce  the  nonlinear  analyses  required  for  plastic 
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yielding  of  aluminum.   To  provide  increased  accuracy, 
each  panel  was  modeled  using  two  meshes.   The  results 
obtained  for  the  coarse  and  fine  meshes  were  extrapolated 
to  a  final  result  using  the  Richardson  extrapolation 
technique  [Refs.  5  and  6]. 

Along  with  ADINA,  a  preprocessor  program,  PSAP1  [Ref .  7] , 
was  used  to  verify  mesh  connectivity  prior  to  analysis  by 
ADINA.   PSAP1  provides  a  graphical  output  of  the  finite 
element  mesh  and  was  coded  for  the  CALCOMP  plotter  installed 
at  NPS  prior  to  1978.   For  this  thesis  PSAP1  was  adapted 
for  use  with  the  VERSATEC  plotter  now  installed  at  NPS. 

The  stress-strain  material  properties  of  the  7075-T6 
aluminum  actually  used  to  make  the  panels  had  to  be  estab- 
lished to  provide  an  accurate  material  model  for  use  with 
ADINA.   Material  testing  was  conducted  to  establish  the 
Young's  Modulus  (E) ,  Poisson's  Ratio  (v) ,  yield  stress  (a  ), 
strain  hardening  modulus  (E  )  and  the  Ramberg-Osgood  co- 
efficients ;3  and  n. 

Comparisons  were  made  to  other  works,  in  addition  to 
the  experiments  conducted  at  NPS.   The  initial  analysis 
involved  a  comparison  of  FEA  to  the  results  of  Howland 
[Ref.  8],  for  a  circular  hole  in  a  finite  strip,  to  vali- 
date the  methods  used.   A  comparison  of  FEA  to  plane 
stress,  slip-line  theory,  for  rigid-perfectly-plastic 
material  was  also  included  as  a  validation  for  the  plastic 
analyses. 
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II.   MATERIAL  PROPERTY  TESTING  OF  7075-T6  ALUMINUM 

The  elastic  and  plastic  material  properties  of  the 
aluminum  panels  were  established  by  tensile  tests  of  uni- 
axial specimens  made  from  the  same  mill  run.   The  speci- 
mens  were  manufactured  and  tested  according  to  current 
ASTM  standards  [Ref.  9].   MICRO-MEASUREMENTS,  EA-13-125AD- 
120,  precision  strain  gages  with  a  temperature  compensated 
bridge  circuit  were  used  on  all  specimens.   Transverse 
gage  sensitivity  errors  were  corrected  according  to  the 
manufacturer's  recommendations  [Ref.  10].   Critical  cross- 
section  measurements  were  made  with  a  micrometer. 

A.   TESTS  FOR  AXIAL  LOADING 

Initial  tests  of  the  two-gaged  specimen,  Fig.  1,  in 
the  MTS  testing  machine  indicated  a  significant  bending 
moment  was  being  produced  by  the  30,000  lb  GRIFF  grips. 
To  investigate  this  problem  further,  tests  were  conducted 
on  both  the  MTS  and  RIEHLE  test  machines  with  a  five- 
gaged  specimen  shown  in  Fig.  2.   The  results  of  these 
axial  loading  tests,  shown  in  Table  I,  verified  that  the 
GRIFF  grips  on  the  MTS  test  machine  do  not  give  axial 
loading.   An  inspection  of  the  gripped  region  on  the  speci- 
men indicated  that  the  jaws  of  the  grip  were  not  applying 
a  uniformly  distributed  force  and  thereby  induced  a  bending 
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moment  by  off -axis  loading  as  shown  in  Fig.  2.   The  grips 
on  the  RIEHLE  test  machine  gripped  evenly  and  a  uniform 
strain  distribution  resulted  as  seen  in  Table  I. 

B.   CHARACTERISTICS  OF  7075-T6  ALUMINUM  PANELS 

The  following  characteristic  properties  were  determined 
from  the  four  specimens  tested. 

1.  Young's  Modulus  (E) 

Tests  were  conducted  using  the  specimen  shown  in 
Fig.  3  on  the  MTS  test  machine  with  10,000  lb  INSTRON  grips, 
which  gripped  the  specimen  evenly.   The  results  of  testing 
three  specimens  are  shown  in  Tables  II  to  IV.   Linear 
regression  in  the  elastic  range  of  all  the  test  results 
determined  a  Young's  Modulus  of  10.12  x  10   psi,  with  a 
correlation  coefficient  of  0.9996. 

2 .  Poisson's  Ratio  (v) 

Tests  were  conducted  using  the  specimen  shown  in 
Fig.  1  on  the  RIEHLE  test  machine  with  10,000  lb  RIEHLE 
grips.   The  results  are  tabulated  in  Table  V.   Linear 
regression  of  these  results  in  the  elastic  region  deter- 
mined Poisson's  Ratio  to  be  0.3256  with  a  correlation 
coefficient  of  0.99996. 

3 .  Yield  Stress  and  Strain  Hardening  Modulus 
These  values,  required  for  ADINA's  bi-linear 

material  model,  were  determined  graphically  using  the 
data  from  Tables  III  and  IV.   Plastic  region  data  in 
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Table  II  is  not  reliable  because  of  excessive  creep  en- 
countered during  that  test. 

0.2%  offset  yield  stress,  a      =  76,000  psi 
strain  hardening  modulus,  E   =  5  66,000  psi 
The  graphical  fit  of  these  values  to  the  test  data  can  be 
seen  in  Fig.  4. 

4 .   Ramberg-Osgood  Coefficients 

The  Ramberg-Osgood  equation  for  elastic-plastic 
stress-strain  characterization  is  given  by: 

c  =  |+6(|)n  (1) 

where : 

£  =  strain 

a  =  stress 

E  =  Young's  modulus. 
The  3  and  n  coefficients  were  determined  graphically  from 
the  data  of  Table  IV,  by  the  method  given  by  Rivello  [Ref.  11] 
The  data  in  Table  IV  gave  the  following  values  which  are 
the  best  fit  to  the  combined  test  data 

S  =  1.479  x  1043 

n  =  21.58 
The  graphical  fit  of  these  values,  in  Eq.  (1) ,  with  the 
test  data  is  also  shown  in  Fig.  4. 
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III.   MODIFICATION  TO  GRAPHICAL  PREPROCESSOR 

The  use  of  a  graphical  preprocessor  program,  such  as 
PSAP1,  is  vital  in  detecting  mesh  errors  that  may  other- 
wise go  unnoticed.   Establishing  the  proper  node  locations 
and  element  geometry  prior  to  analysis  for  a  complex  code 
such  as  ADINA  is  of  utmost  importance. 

A.   PSAP1  MODIFICATIONS 

The  program  PSAP1  was  originally  coded  in  FORTRAN  IV 
for  use  on  the  NPS  IBM  360/370  installation  with  the 
CALCOMP  Model  765  drum  plotter.   The  CALCOMP  system  once 
installed  at  NPS  used  the  +Y  axis  as  the  unlimited  plotting 
axis,  see  Fig.  5.   The  entire  plotting  logic  in  PSAP1  uses 
this  orientation  of  axes  to  allow  multiple  plots  in  a  con- 
tinuous strip.   With  the  VERSATEC  Model  8222A  electrostatic 
plotter  now  installed  at  NPS,  the  +X  axis  becomes  the  un- 
limited plotting  axis,  shown  in  Fig.  5.   To  avoid  an  ex- 
tensive recoding  of  PSAP1  for  use  with  the  VERSATEC  plotter, 
a  simple  coordinate  transformation  of  the  plot  was  made  in 
a  limited  number  of  short  subroutines.   First,  all  installa- 
tion dependent  plotting  calls  used  in  PSAP1  were  identified. 
These  involved  seven  plotter  functions  for  which  new  sub- 
routines were  coded. 
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Function  New  Subroutine 

Initialize  Plotter  CALCMP 

Move  Plotter  Pen  CALPLT 

Letter  on  Plot  NOTATE 

Number  on  Plot  CALNUM 

Determine  Current  Pen  Location  CALWH 

Draw  a  Line  CALINE 

Stop  Plotter  PSTOP 

The  subroutines  listed  above  merely  rotates  the  plot 
to  coincide  with  the  VERSATEC  axis  orientation  and  retain 
all  features  originally  in  PSAP1.   Since  all  plotter 
hardware  code  is  now  isolated  in  these  seven  subroutines, 
future  adaptations  to  other  plotting  systems  is  simplified. 
To  provide  documentation  of  this  update  to  PSAP1,  a  complete 
listing  of  the  new  program  is  provided  in  Appendix  D. 

B.   USE  OF  PSAP1 

Previous  use  of  PSAP1  on  the  IBM  360  system  necessitated 
use  of  a  load  module  since  PSAP1  took  over  one  minute  of 
CPU  time  to  compile.   With  the  new  IBM  3033  system 
compilation  requires  eight  seconds;  however,  use  of  a  load 
module  or  disk  stored  source  code  is  still  recommended 
since  PSAP1  contains  roughly  2,500  lines  of  code.   Appendix 
A  contains  sample  JCL  to  use  PSAP1  on  the  IBM  303  3  MVS 
system.   With  minimal  effort  PSAP1  could  also  be  set  up 
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for  use  on  the  IBM  VM/370  system.   The  user's  manual  for 
PSAP1  is  in  Ref.  7.   In  addition  to  a  mesh  plot,  PSAP1 
provides  a  listing  of  the  node  coordinates,  element 
connectivity  and  several  key  input  values  used  in  execu- 
tion of  ADINA.   This  information  provides  a  useful  check 
of  the  input  data. 
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IV.   FINITE  ELEMENT  ANALYSIS  (FEA) 

A.   DESCRIPTION  OF  MESHES  USED 

1 .   Length  to  Width  Ratios  and  Boundary  Conditions 
Initial  finite  element  models  of  specimens  had 
length  to  width  ratios  near  one,  as  used  by  Garske  [Ref.  12] 
for  his  FEA,  but  they  did  not  provide  the  desired  uniform 
distribution  at  the  loading  boundary.   Specimen  length  to 
width  ratios  of  3-5  were  used  by  Armen,  Pifko  and  Levine 
[Ref.  13]  in  their  FEA  and  by  Stenstrom  [Ref.  1]  in  his 
photoelastic  experiments.   The  criteria  established  to 
determine  uniform  boundary  stress  distribution  was  uni- 
formity in  nodal  displacements  along  the  loaded  edge  as 
discussed  by  Segerlind  [Ref.  14].   In  the  models  used  for 
FEA  in  this  thesis,  nodal  displacements  were  uniform  to 
within  0.1%,  and  the  resulting  stress  distribution  was 
uniform  axially  to  within  0.1%  at  the  panel  ends.   In  all 
cases  two-dimensional,  eight  noded,  isoparametric  elements 
were  used.   These  higher  order  elements  cannot  be  loaded 
in  an  "intuitive"  manner  as  discussed  by  Zienkiewicz 
[Ref.  15,  p.  223].   Figure  6  shows  the  nodal  loading 
required  to  obtain  a  uniform  surface  load. 
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2 .   Element  Meshes 

Two  meshes  were  developed  for  each  panel  analyzed. 
A  reasonable  effort  was  made  to  keep  element  corner  angles 
as  close  to  90   as  possible  to  reduce  the  effect  of  element 
distortion  discussed  by  Hopkins  and  Gifford  [Ref .  16] . 
All  meshes  modeled  a  quarter  of  the  actual  panel  by  using 
the  two  axes  of  symmetry  as  is  common  practice  in  FEA. 
The  step  from  course  to  fine  element  meshes  was  made  so 
that  each  element  in  the  course  mesh  was  subdivided  into 
four  smaller  elements  of  the  same  type.   Such  a  mesh  sub- 
division can  be  expected  to  give  monotonic  convergence  of 
results,  Cook  [Ref.  17],  and  allow  extrapolation  to  results 
of  an  infinitely  fine  mesh.   Figures  7  through  13  illus- 
trate the  element  meshes  used  in  this  analysis  as  plotted 
by  PS API. 

B.   COMPUTATIONAL  PROCEDURES 
1.   Using  ADINA 

Once  the  mesh  has  been  developed,  input  data  is 
prepared  in  accordance  with  the  ADINA  user's  manual  [Ref.  4] 
This  same  set  of  data  is  then  used  as  input  for  PSAP1  to 
check  for  errors  and  provide  a  graphical  display  of  the 
element  mesh.   After  preprocessing  by  PSAP1,  the  data  is 
entered  into  ADINA  for  analysis.   In  the  case  of  linear 
analysis,  two  types  of  stress  output  may  be  specified, 
nodal  point  or  Gauss  integration  point.   Nodal  point  output 
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can  be  computed  for  up  to  eight  node  point  stresses  for 
each  element.   Since  2x2  Gauss  integration  was  used,  four 
Gauss  point  stresses  were  computed  for  each  element.   The 
2x2  Gauss  integration  is  recognized  as  the  most  efficient 
integration  order  for  this  type  of  analysis  [Ref.  15,  p.  284] 
The  linear  analysis  used  an  isotropic  linear  elastic  material 
model  (MODEL  "1"  in  Ref.  4)  which  required  input  of  E  and 
v  material  properties.   The  nonlinear  analysis  allows  only 
Gauss  point  stress  outputs  and  uses  a  bilinear  elastic- 
plastic  material  model,  with  von  Mises  yield  condition  and 
isotropic  strain  hardening  (MODEL  "8"  in  Ref.  4) . 

For  static  analyses  ADINA  uses  a  time  function 
method  to  apply  loads  in  steps.   Linear  analysis  loading 
was  accomplished  in  a  single  step  to  a  nominal  value  of 
3,000  lbs  load.   Nonlinear  analysis  loads  were  applied  in 
ten  steps  to  a  maximum  value,  matching  the  experimental 
loads,  and  then  unloaded  to  zero  in  ten  steps  to  obtain 
residual  stresses.   The  stress  output  from  ADINA  is  a 
listing  of  nodal  or  Gauss  point  stresses  for  each  element. 
Since  the  only  area  of  interest  in  this  analysis  was  the 
distribution  of  stresses  along  the  reduced  cross-section, 
no  large  post-processing  program  was  developed  or  used. 
All  final  computations  using  ADINA  output  data  were  accom- 
plished on  a  HEWLETT-PACKARD  9830A  calculator,  using  short 
programs  coded  in  BASIC.   If  more  extensive  stress 
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distribution  information  were  desired,  some  form  of  auto- 
mated post-processing  would  be  necessary  to  reduce  the 
computational  workload.   At  a  minimum,  nodal  stress  outputs 
by  ADINA  must  be  averaged  to  obtain  unique  values  of  stress 
at  nodes  shared  by  more  than  one  element. 
2 .   Richardson  Extrapolation 

The  use  of  course  and  fine  meshes  allows  extrapola- 
tion to  an  infinitely  fine  mesh  as  discussed  earlier. 
Richardson  extrapolation  [Ref.  5]  was  used  in  this  analysis 
where: 

W    "  w 

aavf„n  =  — - — ~ — —  (2) 

extrap       .  m   ,  m 
n_   -  h 
F      c 

where 

o  -   extrapolated  solution 

extrap        c 

a  =  solution  obtained  with  h=h 

c  c 

a  =  solution  obtained  with  h=h 

r  r 

h       =  linear  dimension  of  course  element 
c 

h„      =  linear  dimension  of  fine  element 

r 

m       =2  (for  this  analysis) 
The  exponent  m  is  determined  by  the  order  of  the  discre- 
tization error  0(h).   Since  h  represents  the  length  of  an 

2 
element  the  element  area  is  represented  h  .   In  a  two 

dimensional  problem  such  as  this  0(h)  is  of  the  order  of 

2 

h  ,  the  area  of  an  element.   In  the  mesh  refinement  scheme 
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1       hF    1 

used  h_  =  -s-h    or  r—  =  ■*■  .  Equation  (2)  can  be  rewritten 

F    2  c    h    2  ^ 
c 

hF  2  hc  2 

*c  <]£  "  *F  <IT> 

0     4.          =   ^U 1~ (3) 

extrap        h  2    hr  2 

c  c 


thus 


aF  "  I  ac 
a  .    =  — — ^ — -  (4) 

extrap       3 

4 


Equation  (4)  then  becomes  the  relation  to  obtain  extrapolated 
stresses  from  coarse  and  fine  mesh  results  in  a  two  dimen- 
sional analysis.   Better  extrapolations  can  be  obtained  by 
using  three  or  more  refined  meshes,  but  the  computational 
effort  increases  significantly. 

3.   Optimal  Stress  Locations  and  Local  Smoothing 

It  is  generally  accepted  that  the  most  accurate 
sampling  points  for  stresses  are  the  Gauss  integration 
points  within  the  element  [Ref.  15,  p.  281,  and  Ref  18]. 
In  this  analysis,  the  nodal  points  are  of  the  greatest 
interest;  thus  a  technique  of  local  smoothing  must  be 
applied  to  the  integration  point  stresses  to  obtain  nodal 
stresses  as  reported  by  Hinton  and  Campbell  [Ref.  19] . 
The  formulation  of  this  local  smoothing  technique  for  ADINA 
elements  is  developed  in  Appendix  B.   The  nodal  values 
obtained  must  then  be  averaged  if  shared  by  two  or  more 
elements. 
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4 .   Computational  Times 

Because  of  the  extremely  large  size  of  ADINA  (about 
17,000  lines  of  code  in  the  NPS  version)  and  the  out  of  core 
solver,  it  does  not  adapt  well  to  time  sharing  systems. 
Using  the  IBM  360  system  at  NPS,  ADINA  required  31  user 
defined  overlays  to  create  a  manageable  load  module  in 
about  30  minutes  of  CPU  time.   With  the  new  IBM  3033  MVS 
system  at  NPS,  ADINA  is  compiled  without  overlays  in  about 
one  minute.   When  using  a  load  module,  the  program  execu- 
tion took  less  than  2  minutes  CPU  time  on  the  IBM  3033. 

In  addition  to  ADINA,  the  preprocessor  (PSAP1) 
and  post-processing  techniques  involve  considerable  time 
and  effort.   Figure  14  is  a  flow  chart  of  the  computational 
procedure  used  in  this  analysis.   An  example  of  the  JCL 
to  use  ADINA  in  load  module  form  on  the  IBM  3033  MVS 
system  v/ith  use  of  the  mass  storage  facilities  is  shown 
in  Appendix  C. 
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V.   RESULTS  OF  ANALYSIS 

A.  CIRCULAR  HOLES  IN  LINEAR  MATERIAL 

The  FEA  results  for  a  circular  hole  in  a  finite  width 
strip  were  used  to  validate  the  elastic  computational  pro- 
cedure discussed  earlier.   The  results  of  Howland  [Ref.  8] 
were  compared  to  both  the  Gauss  point  smoothed  results  and 
the  nodal  output  results  in  Fig.  15.   The  stress  concentra- 
tion factor  a/a      is  referenced  to  the  far-field  stress. 
The  smoothed  results  give  the  best  match  to  the  results 
of  Howland  at  the  edge  of  the  hole,  and  the  only  signifi- 
cant variation  between  the  two  FEA  methods  occurs  within 
the  first  0.25  inches  from  the  edge.   In  order  to  obtain 
the  0.25  inch  stress  value  for  the  coarse  mesh,  in  the 
Gauss  point  smoothed  result,  a  midside  node  value  had  to 
be  obtained  by  the  averaging  method  discussed  in  Appendix 

B.  The  linear  distribution  of  smoothed  stresses  along  the 
sides  of  the  element,  [Ref.  19] ,  appears  to  produce  a 
less  accurate  result  in  this  area  of  extreme  stress 
gradient,  when  compared  to  ADINA's  nodal  output  result. 
This  tendency  was  noted  in  all  cases;  however,  the  peak 
stress  values  from  smoothed  results  consistently  gave 
better  correlation  with  other  investigators  [Ref.  20]. 

A  circular  hole  with  A=0.2  5  was  also  analyzed  and 
compared  to  the  experimental  results   of  Stenstrom  [Ref.  1] 
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along  with  an  interpolation  of  Howland's  results.   The  oa 
experimental  data  correlates  well  with  the  FEA  results; 
however,  the  a   experimental  data  shows  significant  varia- 
tion between  0.125  and  0.375  inches  from  the  edge  of  the 
hole,  as  seen  in  Fig.  16. 

B.   OPPOSITE  U  NOTCHES  IN  LINEAR  MATERIAL 

The  results  for  linear  analysis  are  presented  in  non- 
dimensional  stress  concentration  form;  however,  the 
normalizing  stress  changes.   For  U  notches  KT  is  the 
stress  concentration  factor  referenced  to  the  theoretical, 

nominal  stress  (a  )  in  the  reduced  cross-section  where 

n 

a  =  Load/Area  of  Reduced  Cross-Section. 
1.   Shallow  Notch  Panel 

The  FEA  results  plotted  with  the  experimental 
data  of  Stenstrom  are  shown  in  Fig.  17.   Once  again  both 
FEA  results  are  shown  and  the  variation  for  the  two  methods 
occur  within  0.25  inches  from  the  notch  edge;  however, 
there  was  less  variation  than  was  seen  in  the  circular 
hole  analyses. 

For  this  panel  the  experimental  data  appear  to  be 
uniformly  below  the  FEA  results  for  aft.   The  a   data  shows 
significant  variation  at  the  0.625  inch  point  but  follows 
the  proper  trend  within  0.5  inches  from  the  notch  edge. 
According  to  data  collected  by  Peterson  [Ref .  20] ,  this 
notch  geometry  should  yield  a  maximum  K   =  2.74.   The 
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Gauss  point  smoothed  results  matched  this  value  exactly. 
Stuart  [Ref.  3]  reported  a  K   =  2.6  9  for  this  same  notch 
geometry,  with  a  standard  deviation  of  0.187  in  the  14 
samples  he  measured  by  photoelastic  methods.   Early  photo- 
elastic  work  by  Frocht  [Ref.  21]  determined  K   =  2.7  for 
this  notch  geometry.   The  FEA  results  appear  to  be  in  good 
agreement  with  other  investigators,  for  this  notch  geometry, 
2 .   Deep  Notch  Panel 

Results  of  the  FEA  for  a  deep  notch  were  plotted 
with  Stenstrom's  experimental  data  in  Fig.  18.   The  results 
of  the  two  FEA  methods  again  diverge  within  0.5  inches  from 
the  notch  edge.   FEA  stress  values  at  the  0.25  inch  point 
have  spread  farther  apart  in  this  case  since  the  stress 
gradient  is  very  severe  at  that  point.   The  experimental 

data  correlates  well  for  both  a„  and  o    ;  however,  the  maxi- 

9      r 

mum  experimental  cu  is  considerably  lower  with  a  K   =  3.83. 
Results  reported  by  Stuart  for  this  notch  geometry  was  a 
K  =  4.05  with  a  standard  deviation  of  0. 219  for  14  speci- 
mens measured  by  photoelastic  methods.   Frocht  [Ref.  22] 
reported  a  photoelastic  K =  3.9  for  this  notch  geometry, 
but  concluded  that  the  result  was  5-10%  low,  giving  a 
corrected  range  of  K   from  4.1  to  4.3.   The  Gauss  smoothed 
FEA  result  gave  a  K   =  4.24  which  compares  well  to  an 
empirical  relation  given  by  Peterson  [Ref.  20]  for 
r/d  <  0.25. 
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K^CL-—)  (0.78  +  2.243/t7F) 


2t  2t  2  2t  3 

0.993  +  0.18  ^  -  1.06(^)      +  1.71(^)        (5) 


where 

t   =   notch   depth    (3.9375) 

r  =  notch  radius  (0.625) 

d  =  minimum  width  (15.625) 

D  =  maximum  width  (23.5) 
Inserting  the  values  above  into  Eq.  (5)  gives  a  K   =  4.26, 
which  is  1/2%  above  the  FEA  result.   Other  FEA  results 
reported  by  Armen,  Pifko  and  Levin  [Ref.  13]  and  Griff is 
[Ref.  23],  using  linear  strain  triangle  (LST)  elements, 
produced  K   values  within  5%  of  those  produced  by  use  of 
Eq.  (5) .   The  rectangular  elements  used  in  this  analysis 
are  known  to  give  better  results  than  LST  elements  as 
noted  by  Clough  [Ref.  24].   It  is  clear  that  the  FEA 
results  obtained  for  this  notch  are  in  good  agreement 
with  other  works. 

C.   OPPOSITE  U  NOTCHES  IN  NONLINEAR  MATERIAL 

The  analysis  for  loading  into  the  plastic  region  of  the 
7075-T6  aluminum  was  made  using  the  bilinear  material 
model  discussed  earlier.   The  loads  used  were  selected  to 
match  those  used  in  the  experiments  of  Stenstrom;  thus 
allowing  direct  comparison.   The  strains  obtained  in  those 
experiments  were  used  to  solve  for  stresses  by  use  of  the 
Prandtl-Reuss  plastic  flow  equations. 
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1.   Shallow  Notch  Panel 

The  results  of  FEA  for  the  three  load  cases,  60,000, 
65,000  and  70,000  lbs  are  presented  along  with  the  experi- 
mental results  in  Figs.  19  through  21.   The  afl  results 
compare  well  although  no  trend  for  peak  g~  stress  away  from 
the  notch  edge  is  shown  in  the  experimental  data.   In  all 
cases  the  FEA  determined  the  peak  afl  stress  to  occur  near 
the  yield  boundary,  and  the  gradient  of  the  a„  stress  to 
fall  off  dramatically  in  the  plastic  zone.   This  character- 
istic behavior  of  the  a~  stress  was  reported  by  other 
investigators  [Refs.  13  and  23]  using  FEA  on  2024-T3  aluminum, 
Plane  elastic-plastic  stress  distributions  reported  by 
Frocht  [Ref.  25]  show  similar  trends.   The  experimental 
data  also  shows  a  marked  change  in  the  gradient  of  afi  stress 
within  the  plastic  region.   The  growth  of  this  plastic 
region  is  approximated  using  the  FEA  results  for  this  notch 
in  Figs.  22  through  24. 

Experimental  data  for  the  a   stress  distribution 
matches  the  FEA  results  closely  except  at  the  notch  edge 
where  the  measured  a   does  not  go  to  zero  as  it  should. 
The  characteristic  peak  value  of  a      near  the  plastic 
boundary  as  seen  in  the  FEA  results  is  also  shown  by  the 
test  data. 

The  FEA  residual  stress  computed  upon  unloading 

from  the  three  load  cases  are  shown  in  Figs.  25  through  30. 

The  characteristic  distributions  of  the  a _  residual  stress 

y 
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agrees  with  those  reported  by  others  [Refs.  25,  26  and  27], 
The  experimental  residual  stress  distributions  reported  by 
Stenstrom  show  similar  trends  but  significant  variations 
when  compared  to  the  FEA  results. 
2.   Deep  Notch  Panel 

Three  load  cases  were  computed  to  plastic  loading 
levels;  however,  only  limited  experimental  results  are 
available  for  this  notch  as  seen  in  Fig.  31.   The  30,000  lb 
load  is  just  at  the  onset  of  yield  in  the  notch  root  area. 
Limited  residual  experimental  data  [Ref.  3]  was  available 
for  comparison  in  Figs.  32  and  33  which  are  plots  of  the 

residual  a      and  a       stress  distributions  as  a  result  of  the 

0       r 

three  loading  cases.   The  30,000  lb  load  has  caused  yielding 
in  a  small  region  at  the  root  of  the  notch  as  seen  in  Fig.  34 
Figures  35  and  36  illustrate  that  the  plastic  zone  does  not 
grow  to  the  extent  it  did  in  the  shallow  notch.   The  stress 
gradients  are  very  severe  in  the  deep  notch  and  as  a  result 
the  0.25  inch  sampling  points  used  in  the  FEA  may  be  useful 
in  only  showing  gross  trends  close  to  the  notch  edge.   The 
trends  appear  to  be  much  the  same  as  in  the  shallow  notch, 

with  the  peak  o      and  j   stresses  occurring  near  the  yield 

0       r 

boundary.   The  yield  boundaries  shown  in  Fig.  31  are 
approximations  based  on  qualitative  analysis  of  the  finite 
element  results.   The  <j   experimental  data  for  the  30,000 

0 

lb  load  case  correlates  well  with  the  peak  stress,  again 
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appearing  low  as  it  did  in  the  linear  analysis.   The  residual 

oa    stress  distribution  in  Fig.  32  follow  much  the  same  trends 

as  seen  in  the  shallow  notch  but  with  much  higher  gradients 

within  the  first  0.5  inches  from  the  notch.   The  FEA  results 

for  the  residual  a      stresses  shown  in  Fig.  33  indicate  a 

limitation  in  the  element  size  used  since  the  a   value  of 

r 

the  notch  edge  does  not  return  to  zero  as  it  should.   Because 
of  this  problem  the  data  may  be  questionable  for  showing 
proper  trends  in  the  first  0.5  inch  from  the  notch  edge. 
3 .   Rigid-Perfectly  Plastic  Panel 

A  stress  distribution  for  the  theoretical  material 
used  in  slip-line  theory  was  desired.   By  using  ADINA's 
bilinear  material  model  such  a  material  could  be  approxi- 
mated reasonably  well.   For  this  analysis  a  Young's  Modulus 

2  ft 
(E)  of  10   psi  was  used  to  model  perfect  rigidity.   The 

strain  hardening  modulus  (E  )  was  set  to  zero  to  model 

perfect  plasticity.   Poisson's  ratio  (v)  was  0.4999999, 

2  ft 
as  close  to  0.5  as  the  computer  would  allow.   The  E  =  10 

also  represents  a  computer  limit  in  approaching  an  infinitely 

large  E.   Figure  37  illustrates  the  results  obtained  and 

compares  them  to  a  slip-line  solution.   The  results  are 

normalized  to  the  arbitrary  73,000  psi  yield  stress  used 

in  the  analysis.   The  aQ  values  obtained  agree  exactly 

with  slip-line  theory.   Conversely  the  a   results  do  not 

reflect  the  same  values  as  slip-line  theory,  but  do  show 
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a  similar  trend.   The  growth  of  the  plastic  zone  obtained 
is  shown  in  Figs.  38  through  40. 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  results  obtained  from  the  FEA  have  proven  useful 
in  determining  the  validity  of  experimental  data  gathered 
by  photoelastic   techniques.   The  FEA  results  varied  by 
less  than  1%  when  compared  to  published  analytical  results 
and  handbook  values.   Experimental  data  from  Stenstrom's 
photoelastic  work  correlated  well  with  the  FEA  results. 
The  primary  exception  was  the  residual  stress   experimental 
values,  which  varied  significantly  from  the  FEA  results. 
The  variation  may  be  in  transforming  the  residual  strains 
measured  photoelastically  into  stresses  for  comparison 
with  the  FEA  results,  since  ADINA  only  provides  a  stress 
output.   Limitations  of  ADINA1 s  bilinear  material  model, 
initially  considered  severe,  do  not  appear  to  have  hampered 
this  investigation.   A  possible  exception  is  the  residual 
analyses  where  the  transition  region  from  elastic  to 
plastic  strains  becomes  especially  important.   The  Gauss 
point  smoothed  results  gave  the  best  correlations  at  the 
edge  singularities  in  all  cases;  however,  due  to  the 
limitations  noted  at  0.25  inches  from  the  edge,  the 
results  at  that  specific  point  may  not  be  as  accurate  for 
this  method.   The  nodal  output  results  gave  consistently 
higher  stress  values  at  the  edge  singularities.   Because 
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of  the  severe  stress  gradients  near  the  deep  notch  analyzed 
the  use  of  a  finer  element  mesh  near  the  notch  would  probably 
produce  better  results. 

The  effort  involved  in  developing  two  meshes  such  as 
those  used  in  this  thesis  is  considerable.  An  automatic 
mesh  generation  capability  would  reduce  the  workload  and 
allow  experimentation  with  several  types  of  element  meshes. 

ADINA  proved  to  be  a  useful  and  powerful  program,  as 
expected,  but  something  simpler  and  less  awkward  to  use  may 
be  all  that  is  required  for  two  dimensional  analysis.   Such 
a  system  is  already  in  use  at  NPS  but  does  not  offer  non- 
linear capabilities.   If  use  of  ADINA  is  to  be  continued  in 
this  type  of  investigation,  a  post-processing  program  should 
be  adapted.   There  are  programs  available  to  post-process 
ADINA  data  at  NPS  [Refs.  28  and  29]  but  they  would  require 
modifications  to  work  with  two  dimensional  analyses  and 
the  VERSATEC  plotter. 

Standardized  material  property  testing  would  ease  the 
inevitable  task  of  obtaining  basic  material  properties  for 
use  in  analysis  or  experiments.   Some  form  of  automatic 
data  collection  with  use  of  the  MTS  testing  machine  would 
allow  testing  of  a  larger  sample  population  and  provide 
statistically  more  accurate  information. 
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FIGURE   2 
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CALCOMP    AND   VERSATEC   PLOTTER  AXES 


41 


P=unif  orm     I  oad      over 
element     surface 


nrrm 


i 


2P 
3 

4 


i 


o 


X 


2P 
3 

i 


1 


FIGURE   6 
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FIGURE  8 
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FIGURE    11 
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FIGURE  12 
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FIGURE    19 
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FIGURE  22 
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FIGURE  24 
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FIGURE  27 
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FIGURE  23 
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FIGURE  23 
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FIGURE  31 
DEEP  NOTCH  PLR5TIC  LORDING  RESULTS 
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FIGURE  32 
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F1EURE  33 
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FIGURE  34 


DEEP  NOTCH   30,000  LB  LOAD  PLASTIC  ZONE 
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FIGURE  35 


DEEP  NOTCH   35,000  LB  LOAD  PLASTIC  ZONE 
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FIGURE  36 


DEEP  NOTCH   40,000  LB  LOAD  PLASTIC  ZONE 


72 


m 

u 

cr 

m 


in  .. 


3 
in 

UJ 

cr 


»- 
in 
cc 

.j 

Q. 
I 

>- 


UJ 

u. 

cr 

Ui 
Q. 

I 

Cx 

w 

cr 


•- 

1 

1 

1 1 1 

1 

i 1 1 

>■ 

a 

a 

tr 

a 

Ui 

X 

1- 
u 

2 

a 

\      a 

5ULTS 

-  SLIP  LI 

1 

n 

\       n 

UJ 

cr        ui 

T                        -  • 

uj 

\ 

£C          -J 

UJ 

U.          Ck 

1            -J 

1            g 

a 

\    a 

a 

a 

\  a 

a 

\a 

1 

\ 

a 

i_ 

Na 

\ 

r 

? 

, 

i 

i 

h 

L 
1 

1 

1 

r 

i 

i 1 

i 

1 

C9 
X 


P4 


m 
is 


in 


si 


si 


3 

I 


D/D '  mDUJ  553HJS 


ui 
rn 


23 


ui 


w 


UJ 

z 

s 

m 


UJ 


tr 

a 

z 


UI 


73 


FIGURE  38 


RIGID-PERFECTLY-PLASTIC   INITIAL  PLASTIC  ZONE 
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FIGURE  39 


RIGID-PERFECTLY-PLASTIC  INTERMEDIATE  PLASTIC  ZONE 
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FIGURE  40 


RIGID-PERFECTLY-PLASTIC   FINAL  PLASTIC  ZONE 
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TABLE'  I.        MTS  AND  REIHLE  5  GAGE  TEST  RESULTS 

—6 
All  Strains  are  10    in/in 

MTS  Test  Machine 


LOAD 
lbs 


STRAINS 


1000 

-3 

421 

814 

1235 

1638 

2000 

338 

1003 

1619 

2279 

2903 

3000 

960 

1714 

2431 

3186 

3908 

REIHLE  Test  Machine 


LOAD 

STRAINS 

lbs 

£1 

Z2 

£3 

£4 

£5 

1000 

800 

763 

755 

741 

726 

2000 

1564 

1540 

1546 

1543 

15  35 

3000 

2370 

2334 

2363 

2377 

2388 
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TABLE    II .       MTS    SPECIMEN    A   TEST    RESULTS 


Cross-section  =   0.03975    in' 


Load ,  lbs 


-6 


Strain,  e , 10    in/in 


256 

615 

503 

1204 

750 

1301 

1005 

2423 

1255 

3042 

1505 

3665 

1773 

4352 

2  0  0  3 

4925 

•-■  -j  cr  •- 

5537 

2503 

6355 

C  l  J  -J 

7365 

2905 

8  2  3  0 

2984 

9045 

3037 

1015O 

TABLE    III.       MTS    SPECIMEN   B    TEST    RESULTS 


Cross-section  =   0.0  3975   in 


Load,    lbs 


-6 


Strain, e    ,    10        in/in 


~j  u  y 
759 

100  0 

1250 

1500 

1750 

2  0  0  0 

2250 
2500 
275  0 

2  9  0  0 

3  0  0  0 

3  0  5  0 

3100 

3125 


1368 
2459 
3  6  8  0 

3  7  0  0 
4328 
4958 

5  6  3  0 

6  3  8  0 
7340 
8200 
9  2  5  0 

1  0  2  0  0 

11553 
12300 
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TABLE    IV.  MTS    SPECIMEN    C   TEST    PESULTS 


.     2 


Cross-section  =   0.03975   in 


Load,    lbs. 


—6 
S  train , e , 10      in/in 


,_ 

..  ..  .. 

•-  !  £ 

bb  :; 

503 

1228 

762 

1850 

1  0  0  8 

2451 

1231 

3  9  0  5 

1506 

3639 

1755 

43  I 

2  8  0  0 

4  9  4  0 

ii  J  .J 

5607 

2503 

6  3  2  5 

2750 

.■  '  J  ' ..'  i  j 

i  ._  -...'  U 

2  'f1  0  0 

8128 

3  0  0  0 

9  2  0  0 

3  0  6  0 

i.  0500 

3  1  0  8 

11588 

3125 

!.  2588 

SO 


TABLE    V. 


REIHLE    SPECIMEN    TEST    RESULTS 

2 
Cross-section  =   0.12    in 


LOAD 
lbs. 


STRAIN, ei 
10      in/ in 


STRAIN,£2 

-6 
10      in/in 


5  @  8 

345 

■- 1  i.  5 

1  8  8  8 

728 

-248 

1  5  0  0 

1  1  1  4 

364 

2  8  8  8 

1581 

i-92 

2  5  8  8 

1395 

-614 

3  8  8  8 

2291 

-  , '  4  !~i 

3  5  8  8 

2785 

37  1 

4  8  8  8 

6  y  y  5 

....  Q  Q ;!."] 

4508 

3515 

i.  :.  ..:.  b 

5  8  8  8 

3912 

....  |    :  •;  1=] 

5  5  8  8 

4332 

■•- :.  .  t 

5888 

4758 

- 1  5 :  8 

6588 

■J  i  o  o 

.1.  • .' '  ■•'  , 

7  8  8  8 

■J  J  ?  •-' 

-17S4 

7  5  8  8 

b  8  7  5 

-1923 

3  8  8  8 

6  6  4  9 

-2183 

3  5  8  8 

f  o 8  5 

-  2  3  7'  5 

9888 

8  6  6  3 

-  2  8  8  8 

9588 

12245 

-44  35 
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TABLE   VI. 


X    =    0.2      HOWLAND   DATA 


DISTANCE 
FROM  HOLE 
in. 


9/  a 


0  .  0 

3 .  1  4 

0.  5 

1  .  57 

1.0 

1  .  26 

1.5 

1  ..  L  S 

2  a  0 

1 .  1  1 

2.5 

1  .  07 

3 .  0 

1  .  05 

-   c 

1  .  0  1 

4 .  0 

y.  97 

TABLE   VII 


X   =    0.25    HOWLAND   DATA 


DISTANCE 
FROM  HOLE 
in. 


9/  a 


0  .  0  0  0 

-' '  ;-  '■-' 

0.422 

1.  75 

0 .  3  2  S 

1.30 

1.234 

1.22 

1  ,.641 

L .  .1 3 

2.  047 

1.  U7 

2.453 

I  .04 

2.359 

0.  97 

3.  063 

0.  95 
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TABLE  VIII. 


DISTANCE 


X    =  0.2  FEA  RESULTS 
a     =  5000  psi 


-  NODAL  OUTPUT 


FROM  HOLE,    in. 

Sg/      PSI 

<?r   ,   psi 

0  .  0  0 

15938.2 

;'  5  1  ,  7 

M         •"■'  C 

1Q259.7 

28  75.,  9 

0  .  50 

"7  O  ~>  cr 
i    O  I    ■_'  ,i    i 

1.733.5 

0.75 

6  8  7  7  .  0 

1673,,  1 

1  .  00 

6  3  0  7  „  4 

1334.7 

1.25 

5995. 4 

1037,,  2 

1.50 

5797.7 

3  74,,  4 

1.75 

5633. 9 

7  0  0  .  5 

2  .  0  0 

5509. 8 

5  5  6 .  3 

TABLE    IX. 


X     =0.2    FEA   RESULTS    -    GAUSS    OUTPUT 
aoo    =    5000   psi 


DISTANCE 


FROM  HOLE,    in. 

GQ     'PS1 

a      ,psi 

r 

0. 

0  0 

15595 

4 

L28.  9 

0. 

25 

9432. 

6 

2134.  1 

0. 

50 

7953. 

'•-' 

1826.6 

0. 

-?cr 
<    J 

6  7  8  3 . 

i 

1658. 4 

1  . 

8  0 

•6316. 

1 

1348. 7 

1. 

•-,  cr 
cl  -J 

5975. 

6 

1983.3 

1. 

50 

5795. 

•-• 

2^}„  1 

1. 

5623. 

8 

695.  7 

W.      1 

00 

5499. 

9 

554.  9 

TABLE    X. 

DISTANCE 
FROM  HOLE,    in, 


X    =   0.25    FEA   RESULTS    -  NODAL  OUTPUT 
a     =  4500  psi 
aQ    ,psi  ar   ,psi 


0  .  0  0 

14745. 

.j 

3  2  9  .  S 

0.  25 

9447. 

0 

2131.1 

0  .  5  0 

7232. 

6 

1599. 3 

0.  75 

6  3  0  4  . 

< 

1473. 5 

1  .  00 

5771. 

r 

1139.6 

1.25 

5471. 

2 

391.6 

1.50 

5266. 

9 

671  .  0 

1.75 

5  9  9 1  . 

( 

+  39.  1 

2.  00 

4940. 

- 

354.  5 

TABLE    XI. 

DISTANCE 
FROM  HOLE,    in. 


X    =    0.25    FEA   RESULTS    -    GAUSS   OUTPUT 
a     =  4500  psi 


y . 

uu 

0. 

■-•  cr 

£  .J 

0. 

50 

0. 

-?  cr 
i  J 

1  , 

00 

1. 

■-.  cr 

i—  ■- 

1. 

50 

1. 

-?  cr 
r  J 

2  .  0  0 


,psi 


14  422. 3 
3721.2 
7317.  1 
6213.6 
5779. 6 
5453. 3 
5267.  1 
5985.9 
4937. 9 


a   ,psi 

r 


927.  1 


1452.4 

1151.5 

3  9  1 .  8 

436.9 


3  4 


TABLE  XII.    SHALLOW  NOTCH  FEA  LINEAR  RESULTS  -  NODAL 

=  4225  psi 


n 


DISTANCE 
FROM  NOTCH,    in. 


-psi 


a   ,psi 
r  r 


0  .  6  0 

11534.9 

94,5 

S..25 

896"6.  1 

1279,9 

0.  50 

7335. 6 

1501 , 9 

0.  75 

6347. 3 

.1731.7 

1  .  00 

5655. 0 

1  7 1  3 .  2 

1.25 

5134,  1 

1703.5 

1  .  59 

4832.4 

1626. 7 

1.75 

4  564.7 

1 568, 6 

2 .  0  8 

4361,  1 

1434, 7 

TABLE  XIII.   SHALLOW  NOTCH  FEA  LINEAR  RESULTS  -  GAUSS 


a   =  4225  Dsi 
n 


DISTANCE 


FROM  NOTCH,    in. 

CTg      /PSi 

ar   /Psi 

8. 

88 

11538.2 

11.5 

8 . 

"  c 

3727.  1 

1177.7 

8 . 

58 

"7  •-!  er  cr      cr 

i    .'<  .j  .j  ,  .j 

1594,3 

0  . 

1      J 

6  2  3  8 .  9 

1713.6 

1. 

8  8 

5661  . 2 

[ 72Q. 4 

1. 

•-■  cr 

5159. 3 

1699. 7 

1. 

58 

4334.6 

1634. 2 

1. 

i    J 

4552.3 

1561.7 

—  • 

88 

4353.6 

1436. 0 
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TABLE  XIV. 


DEEP  NOTCH  FEA  LINEAR  RESULTS  -  NODAL 


DISTANCE 
FROM  NOTCH,  in. 


a     =   4800  psi 
n 


,psi 


<?r     /PSI 


0  .  0  0 

21  142.3 

-:  cr  c   -i 

i  J  -J  „  ii 

0/25 

1 1 929 „ 9 

-t  4  O  '-!•  ,.  O 

0 .  5  0 

3436,,  3 

353U.4 

8.  75 

6  9  7  6 .  6 

33  3- 1.  4 

1  .  30 

6876. 1 

c>  i  -J  4  .  4 

1.25 

5527. 0 

2383.4 

1.50 

5139.2 

2623.4 

1.75 

4342. 4 

24  12,,  0 

2  .  0  0 

4617.3 

2234. 2 

TABLE    XV. 


DEEP   NOTCH    FEA   LINEAR   RESULTS    -    GAUSS 


DISTANCE 
FROM  NOTCH,    in, 


a      =   4800   psi 

n  e 


,psi 


or    rPSX 


0  .  0  0 

2  0  3  6  3  . 

1 

314.  0 

0.25 

10532. 

.j 

h  j".  ■",  vsr       -, 
>-t  ■■_]  r~  j  „  .J. 

0  .  5  0 

3  6  0  6  . 

0 

3611. 3 

0.75 

6352. 

4 

3501  „ 9 

1  .  00 

6036. 

r  , 

;•::;  1  j  -  (  ■;- 

1.25 

5496. 

1 

2334.7 

1 .  50 

5141. 

,-;, 

264  3 .  iJ 

1.75 

4330. 

1 

2421.2 

2  .  0  0 

4622. 

0 

.-, .-,  -, .:,  ,:, 

86 


TABLE    XVI.       SHALLOW   NOTCH   FEA  NONLINEAR      60,000    lb    LOAD 


DISTANCE 


NO    LOAD    RESIDUALS 


FROM  NOTCH,. 

m.  ag  ,psi 

a   ,psi 

r 

Cg  rpsi 

a   ,psi 

r 

0  .  0  0 

30024. 7 

1994. 5 

-33411.4 

439.,  1 

0.25 

•31333.7 

9523. 3 

-  4  2  5  0  „  0 

-  2  6  9  9  .  1 

0  .  5  0 

31333.6 

16052. 6 

9  0  8  6  „  3 

2133.3 

9 .  75 

35543. 3 

13993.2 

2  9  '4  6 ,  6 

2133.4 

1  .  00 

53431.4 

1 S 799.6 

1  9  9  2 ..  3 

1796. 7 

1.25 

52670.  1 

13347.9 

1167,,  3 

1501  . 6 

1 .  50 

49111.  1 

17477. 1 

3  ui  5  „  ,  ■ 

1 2  5  0  .  8 

1.75 

46030. 7 

16567.3 

394..  4 

1046. 7 

2  .  0  0 

43373.3 

15667. S 

471.  6 

335.  0 

TABLE    XVII.         SHALLOW  NOTCH   FEA  NONLINEAR     65,000    lb      LOAD 


NO   LOAD    RESIDUALS 


DISTANCE 
FROM  NOTCH, in.         a, 


,psi 


ar   ,psx 


rPSl 


a      ,psi 
r 


0.  30 

3  0  8  9  0  .  4 

10O4. 6 

-42992.  : 

224.  7 

0.  25 

3  2  9  3  9 .  3 

10229, 5 

-9796. 3 

-  4  7  8  4  .  y 

0  .  5  0 

8  3  6  4  2 .  7 

16054. 6 

6944,,  3 

-711.2 

0.75 

74622.0 

21475.5 

7507,,  7 

3214.2 

1  .  00 

64770. 4 

21171.3 

:'  3  2  0  .  9 

2 ;-::":  ft.  7 

1.25 

57391.2 

20593. 6 

2202. 3 

2419. y 

1.50 

5  3 766. 7 

19537.3 

1524. 5 

2013.5 

1.75 

50321.  1 

13474.3 

1039.5 

1710.1 

2  .  00 

47943.2 

17421.4 

i~,  tr  ■-.   -? 
O  J  Z<  .  i 

1443.0 

37 


TABLE    XVIII.         SHALLOW   NOTCH    FEA      NONLINEAR      70,0001b      LOAD 


NO    LOAD    RESIDUALS 


DISTANCE 

FROM  NOTCH, in. 


Cg        fPSl 


a      ,psi 
r 


tfg    rpsi 


a      ,psi 
r     e 


0  .  0  0 

31957. 9 

704.5 

-  5  U  6  7  0 .  6 

63.  9 

0 .  2  5 

32322. 6 

1 0  6  6  8 .  8 

-16176.1:; 

-3235. 1 

0  .  '-•  0 

34007. 3 

15514.4 

-376. 9 

-2366. 8 

0 .  75 

37017. 4 

22662. 9 

15132.6 

2762. 3 

1  .  00 

7 1  2: 1  8 .  0 

23765. 2 

5631.  1 

3983. 0 

1.25 

63732. 9 

23215.2 

3763.  1 

3 6 16.6 

1,50 

53941 . 5 

21920. 1 

2  6  6  9  „  9 

3  0  3  0  .  0 

1.75 

54337.3 

20670. 6 

1365,,  4 

2596.7 

2.  00 

52177.0 

19426. 6 

1443.2 

2207. 0 

TABLE    XIX. 


DEEP   NOTCH      FEA      NONLINEAR      30,000    lb      LOAD 


NO    LOAD    RESIDUALS 


DISTANCE 


FROM  NOTCH,  in.  a       ,psi 

9 

a   ,psi 

r 

a9  ,psi 

a   /psi 

r 

0.08 

33521 . 0 

23  32. 4 

-14133. 3 

_  -7  •-,  cr   -, 

0=25 

F;  8@  :-::=:.  6 

20971 . 3 

5630.  1 

1402. 5 

0 ,  5  0 

43309.  1 

13640. 5 

O  O  i~  <  7 

324.  6 

3.75 

34402. 6 

17339.  1 

193.3 

463.9 

1  „  90 

30522.9 

16076. 1 

125.3 

300.5 

1.25 

27543. 7 

14574.3 

33.,  9. 

215.2: 

1.50 

•-  cr  -?  c  — ■   ,-j 

13323.7 

6  7 .  3 

160.  7 

1.75 

24133.3 

12195.0 

53.  5 

126.3 

2  .  0  0 

23141.6 

11265.7 

44.9 

130.9 
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TABLE    XX. 


DEEP   NOTCH      FEA      NONLINEAR      35,000    lb      LOAD 


NO    LOAD    PESIDUALS 


DISTANCE 
FROM  NOTCH, in. 


ae  '?sl 


a      ,psi 

r 


°a     'PS1 


rpsx 


0  .  0  6 

79020. 4 

1253. 

1 

-  :;;  9  fi  6  ■'-'■ .  2 

-19 0  9 . 2 

0 .  2  5 

78215.7 

24593. 

9 

16  343. 7 

1113.1 

0 .  5  0 

5  0  0  9  4 .  6 

d  d  fc  i'"  • j  . 

1 

105.7 

1699. 5 

0,75 

4051b. 0 

21359. 

c 
.J 

624.  1 

1002. 3 

1  .  00 

■Z<  -J  O  4  •-■  .  1 

19156. 

9 

3  8  0  .  6 

716.7 

1.25 

32290. 3 

17299. 

1 

2  ■_'  r  .  3 

521.7 

1.58 

30172. 9 

15732. 

1 

2  0  3  „  5 

4  0  5  .  9 

1.75 

23314.3 

144  10. 

'-' 

156.9 

3 19.5 

2  .  0  0 

27079. 2 

13294. 

■4 

1  3  3  „  r' 

261 .  5 

TABLE    XXI. 


DEEP   NOTCH      FEA      NONLINEAR     40,000    lb      LOAD 


NO    LOAD    RESIDUALS 


DISTANCE 
FROM  NOTCH, in. 


rPSl 


a      ,psi 
r 


,psi 


a      ,psi 
r 


0  .  00 

8 0929. 5 

335.4 

-  5  2  9  9  0  .  9 

-2121.6 

0.25 

8 9730.7 

24933.3 

2  0  5  0  3  „  1 

-  '2 1  -1 2  .  3 

0 .  5  8 

56572.6 

26315.4 

-613, 4 

2977. 3 

0.  75 

46984.6 

248  73. 1 

1  4  6  3  .  2 

1643.  1 

1.00 

41400. 6 

22236. 3 

6  1  0  .  0 

1216.6 

1 .  25 

37164.5 

20063. 3 

555.  9 

392.  7 

1 .  50 

34690. 9 

13237.2 

438.5 

710.5 

1.75 

32512.6 

16672.2 

330.  7 

564.  3 

2.00 

31073. 3 

15363.3 

232.  0 

470.  3 
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TABLE    XXII. 


RIGID   -    PERFECTLY    -    PLASTIC    RESULTS 


DISTANCE 
FROM  NOTCH, in. 


a     =    73,000  psi 

y 


°9     rPSl 


a      ,psi 

r 


0  .  0  0 

73334. 4 

5 1  9  .  9 

0.25 

31730.5 

20863 „ 6 

0  .  5  0 

8  '3252.  6 

293  13,.  3 

0.  75 

34074. 1 

37321. 1 

1  .  00 

34103.  1 

'4  S  2  4  4' .  id 

1   '"' c" 

34091 . 3 

42244,,  3 

1.50 

34051  .  5 

4324  9,.;? 

■i    —•  ET 
i  .  !'  J 

3: 4  0  8  5  .  0 

4  4  y  3  2  .  9 

2.  90 

33353. 1 

■4  3  o  i7  5  .  6 

TABLE    XXIII.         EXPERIMENTAL    DATA    \   =    0.25    HOLE      LINEAR  LOADING 

a     =   10 ,749  psi 


DISTANCE 
FROM  HOLE,    in, 


°Q         »PS1 


a      ,psi 
r 


0  .  0  0  0 

35221.0 

49.  5 

0.  125 

243  3  9 . 0 

-1350,.  5 

0.  250 

19309. 5 

6  6  9  .  5 

0.375 

17799. 0 

2311.5 

9 .  5  0  0 

16588.0 

3422. 3 
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TABLE    XXIV. 


EXPERIMENTAL   DATA      SHALLOW   NOTCH      LINEAR  LOADING 


15,000   lb  Load 


DISTANCE 
FROM  HOLE,  in, 


sQ    ,psi 


^r    /PSl 


0  .  00  8 

23311. 

"7 

Z<  o  ,<  cz. 

0.  125 

22295. 

0 

13 92. 6 

0.  253 

1 9  9  9  9  . 

•-"' 

4  3  5  £  „  7 

0.375 

17277. 

2  3  3  7  .  3 

0  ■  5  0  0 

15793. 

f 

2338. 5 

0 .  325 

13379. 

I 

1929. 3 

TABLE    XXV.    EXPERIMENTAL   DATA  SHALLOW  NOTCH      60,000    lb    LOAD 


NO    LOAD    RESIDUALS 


DISTANCE 
FROM  NOTCH, in. 


,psi 


a      ,psi 
r 


°q     /PSl 


ar    ,psi 


8 .  3  S  0 

3  1 1 1  1  .  5 

55  51  . 8 

-'d'SdA'd  ..  8 

4 .  8 

0.  125 

79710. 7 

3735. 3 

-19337,.  8 

1 7  5  .  0 

0.25  0 

73112.0 

8141.5 

- 1  3  4  8  9  .  8 

-1844. 0 

0.375 

77139.3 

1 1399. 5 

-  1  2  3  2  0  .  3 

-1423. 0 

0  .  5  0  0 

71207. 2 

10539. 4 

-123  1  9 , 0 

-5033  ,  0 

0.325 

71433.4 

192  3  0 . 8 

-  1  8  3  0  3  .  3 

'"j  i-i"  .j  -.j  .  y 

0.  750 

33045.3 

1  3  5  3  0  .  9 

0.375 

59477.4 

17743.3 

1  .  080 

57143.4 

19275.4 
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TABLE    XXVI.  EXPERIMENTAL    DATA   SHALLOW   NOTCH      65,000    lb    LOAD 


NO   LOAD    RESIDUALS 


DISTANCE 
FROM  NOTCH, in. 


aQ    rpsi 


ar   ,psi 


afl    /Psi 


ar    /Psi 


8 .  yyy 

33873. 

Q 

7  0  5  8  .  8 

_  "'  |L"  '"'  ,  1  r~.        |"i 

-  :.  9  1  .  8 

8,  125 

8 8957. 

b 

6  0  7  8  .  9 

-  32 j 4 3  .  U 

-337. 0 

8 .  258 

- '  I™.  I-  cr  i 

i  i  o  •-'  4  . 

9 

7468.  L 

""  2 : ::;  2  0  b  ,,  y 

-1554. 0 

9.  375 

7  9  4  8  2 . 

'-' 

10373. 9 

-  2  3  9  1  4  .  3 

-2381.8 

8  .  5  8  8 

79333. 

4 

166  3 4 . 8 

■•-  2  8  1 9  2 .  8 

-  2  6  8  6 .  0 

8.  325 

r  i  o  O  i'  ■ 

8 

1  9  8  6  7  .  6 

~  1 '-- ':  8  0 .  0 

- 3 3 92.0 

8 .  75  8 

<     C    J  :J   b  . 

4 

197  8 6 >  ? 

8,  375 

34333. 

2 

17433. 1 

1 .  388 

61648. 

3 

19834. 1 

TABLE    XXVII.         EXPERIMENTAL   DATA   SHALLOW   NOTCH      70,000    lb   LOAD 

NO    LOAD    RESIDUALS 
a      ,psi  a 


DISTANCE 
FROM  NOTCH, in 


rpsx 


,psi 


ar    /PSI 


8 .  8  8  8 

34743. 9 

8 569. 8 

-58791 

8 

43.  8 

0.  125 

33492.9 

9873. 9 

-4r  0  :.f:  , 

0 

-736. 0 

8.  250 

34253.3 

13713.2 

-40774. 

8 

-378.0 

8.375 

32346.3 

13815. 8 

-  3  -1-  6  3  8  , 

8 

-191 6 . 0 

8 .  5  8  0 

S  2  5  £  4 .  6 

16314. 9 

-  3  y  8  1  7 . 

8 

-  L  3  1  6  .  0 

0.  625 

30446. 4 

17441  . 6 

-24798., 

0 

-1551 . 0 

0.  750 

80455. 5 

23190. 4 

8.375 

76563. 9 

22923. 7 

1  .  800 

72264. 1 

2  3  604. y 
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TABLE  XXVIII.    EXPERIMENTAL  DATA   DEEP  NOTCH   LINEAR  LOADING 

15,000  lb  Load 


DISTANCE 
FROM  NOTCH, in. 


Oa     'PSl 


a      ,psi 

r 


0 .  S  5  0 

45907. 0 

-45-  u 

0.  125 

35372.0 

7  6  9  3  ..  y 

0.  250 

30410. 0 

11671. 0 

a    ■■■  -  =: 
y .  o  i  _i 

26177. 0 

12395, 0 

0  .  5  0  0 

289 J 9. y 

9  8  7  9  .  y 

0.  625 

13373.0 

9  0  6  y .  y 

0.  750 

16391 . 0 

3934, y 

TABLE   XXIX.  EXPERIMENTAL   DATA   DEEP   NOTCH      30,000    lb    LOAD 

Elastic   -  Plastic 


DISTANCE 
FROM  NOTCH, in. 


rPSl 


a      ,psi 

r 


0  .  0  0  0 

76156. 

6 

•"'  ■=:  '.:•  ■?  a 

0.  125 

6  3  9  9  6 . 

- 

6057, 7 

0 .  25  0 

5  0  047. 

-, 

10456,,  3 

0.  375 

43992. 

■-, 

14435,5 

0  .  5  0  0 

39743. 

ifi 

16435. 6 

8.  625 

3  5  0  3  0 . 

'■-• 

15399. 3 

0.  750 

3  3  2  5  6 . 

;:: 

16725. y 
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APPENDIX  A 
PSAP1   JCL 


STANDARD    JDB    CARD    

//    EXEC    FRTXCLSP 

//FORT.SYSPRINT       DD       DUMMY 

//FORT.SYSIN    DD    UN  IT=3330  »VOL=SER=OI  SK02  t 

//    DSN=S2  939.PSAP(PSAP) , D  ISP=SHR, LABEL=M  » , t  IN) 

//    DD    UNIT=3330tVOL=SE*=OISK02fDSN=S2939.PSAP(PLOT) , 

//    DISP=SHR,LABEL=(  ,,,  IN) 

//    DD    UNIT=3330,VOL=SER=0ISK02,DSN=S2939.PSAP< INIT) , 

//    DISP=SHR,LABEL=( t,,IN) 

//    DD    UNIT=3330,VOL=SER=DISK02,DSN=S2939.PSAP(ELER) , 

//    DISP=SHR,LA3ELM  ,.»IN) 

//    DD    UNIT  =  3330,VOL  =  SE:*  =  DISK02,DSN=S2939.PSAP<SAPF)  , 

//    DISP=SHRTLABEL=( ,,, IN) 

//    DD    UNIT*3330fVOL=SER=DISK02tDSN=S2939.PSAP(ADNA) f 

//    DISP=SHR,LABEL=( ,t» IN) 

//  DD  UNIT=3330,VOL-SER=DISK02fDSN=S2939.PSAP(AUXL) , 

//  DISP=SHR,LABEL  =  (  ,,, IN) 

//  DD  UNIT=3330 ,VOL=SER=DI SK02,DSN= S2939. PSAP< ADPT) , 

//  DISP=SHR,LABEL  =  (  ,,, IN) 

//  DD  * 

C  ***  MAIN  PROGRAM  *** 

DIMENSION    ZZZC3233) ,DI SPD ( 5 ,3 »830  ) 

CALL    PSAPK  ZZZ,32D3,DISPD,300) 

STOP 

END 
*************    DELIMITER    CARD       (/*)     **************** 
//GO.FT10F001    DD    UN IT  =  S YSDA, D  ISP= ( , PASS  )  f 
//    SPACE=(CYL,( 2, 2) ) ,DSN=££TEMP1 

//GO.SYSIM    DD    * 

************* ************************** ****** 

*********  INSERT  PSA'l  DATA  HERE  *********** 
********************************************* 
*************    DELIMITER    CARD       (/*)    **************** 
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APPENDIX  B 

LOCAL  LEAST  SQUARES  SMOOTHING 
I 


•  Corner  Nodes 

o  Midside  Nodes 

+  2x2  Gauss  Points 


Two-Dimensional  Isoparametric  Element  from  ADINA  [Ref .  4] 


The  local  smoothing  expression  from  Hinton  and  Campbell 
[Ref.  19]  in  ADINA  coordinates  becomes 

\ 


(    \ 

c 

B 

B 

A 

32 
S3 

)     = 

B 
A 

A 
B 

C 
B 

B 
C 

94 

B 

C 

A 

B 

^      1 

X 


/3  1 

where  A  =  1  +  — T  ,  B  =  -  ■=-  and  C  =  1  - 


II 


III 


IV 


/3 
2  * 


With  a,,  Qp,  a -.    and  a.  representing  the  smoothed  corner  node 
stresses  and  a  ,  cj__,  aTTT'  anc^  aTV  as  the  unsmoothed 
stresses  at  the  Gauss  integration  points,  this  expression 
can  be  written  in  an  equivalent  form. 


a. 


al/ 


2 

V  ) 


A 

B 

C 

B 

B 

A 

B 

C 

C 

B 

A 

B 

B 

C 

B 

A 

X 


\ 


III 


IV 


II 
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The  midside  node  stress  values  may  be  obtained  by  averaging 
the  values  at  the  associated  corner  nodes,  since  the  distri- 
bution of  the  smoothed  stresses  is  linear  along  the  sides  of 
the  element.   Smoothed  stress  values  obtained  by  this  least 
squares  method  should  subsequently  be  averaged  to  obtain 
unique  values  at  nodal  points  shared  by  adjacent  elements. 
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APPENDIX  C 
ADINA  JCL 


//  EX 
//FOR 
//FOR 


***** 

//IKE 

//LKE 

IN 

EN 
***** 


EC  F 
T.SY 
T.SV 

IMP 
RE  A 
COM 
MAX 
CAL 

sro 

^END 

,v»  -**  1*  -I* 

D.  JS 

D.SY 

CUD 

TRY 

**** 

//G0.FT07 
//  DCB=( 
//GO.FTOl 
//  DC3=( 
//GO.FTD2 
//  OCB=( 
//G0.FTD3 
//  0C8=( 
//G0.FT04 
//  0C3=( 
//G0.FT08 
//  DC3=( 
//GO.FT09 
//  OCB={ 
//GO.FTLO 
//  0C3=( 
//GO.FTL1 
//  DC3=< 
//G0.FT12 
//  DCB=( 
//G0.FTL3 
//  DCB=( 
//G0.FT56 
//  DSN=FO 
//G0.FT57 
//  DC3=( 
//GO.SYSI 

********* 


-  STANDARD 
ORTXCLGtRE 
SPRINT  DO 
SIN  00  * 
LICIT  REAL 
L  A 

MON  A( 1200 
=120000 

L  exe:max 

D 


J33 
0I3N 


CARO    • 
=2000K 
Y 


DJW 

*8( A-H,0-Z) 
10) 

) 


****      [) 

00    0  0 
SIN    00 
E    OS  DO 

'•IAIN 
****    o 

P001  0 
RECFM= 
F001  D 
RECFM= 
F001  D 
RECFM= 
F001  0 
RECFM= 
F001  0 
RECFM= 
F001  0 
RECFM= 
F001  0 
RECFM= 
F001    0 

RECFMr 

F001     D 

RECFM= 

F001    0 

RECFM= 

F001    0 

DECFM= 

F001    0 

099. TE 

«=001    D 

RECFM= 

N  00  * 
****** 


ELlM 

DISP 


ITER 
=  SHR 


CARO       (/*)     **************** 
,DSiN='<1SS.S2939.  AOINA 


(LOAO^ ) 


ELIM 
D  UM 
VBS, 
D  UN 
73  S» 
0  UN 
VBS, 
D  UN 
VBS, 
0  UN 
VBS, 
0  UN 
V3S, 
0  UN 
VBS, 
D  UN 
VBS, 
D  UN 
VBS, 
D  UN 
VBS, 
0  UN 
VBS, 
D  UN 
MP,0 
0  UN 
VBS, 


ITE* 
IT  =  S 
3L<S 
IT  =  S 
BL<3 
IT  =  S 
BL<S 
IT  =  S 
3L<S 
IT  =  S 
BL<S 
IT  =  S 
3L<S 
IT=S 
BL<S 
IT  =  S 
BLKS 
IT  =  S 
3LKS 
IT  =  S 
BL<S 
IT=S 
3L<S 
IT  =  3 
IS°  = 
IT  =  S 
3L<S 


CARO 
YSDA, 
IZE  =  ^ 
YSDA, 
IZE=^ 
VSDA, 
IZE  =  ^ 
YSDA, 
IZE=1 
YSDA, 
IZE  =  1 
/SDA, 
IZE  =  ^ 
YSDA, 
IZE=^ 
YSDA, 
IZE='+ 
/SDA, 
IZE  =  1 
/SDA, 
IZE=1 
YSDA, 
IZE  =  1 
330, V 
SHR,L 
YSDA, 
IZE  =  1 


(/* 

OISP  = 
000)  , 
OISP  = 
000)  , 
DISP  = 
000)  , 
OISP  = 
000)  , 
OISP  = 
000)  , 
DISP  = 
000)  , 
DISP  = 
000)  , 
DISP  = 
000)  , 
OISP  = 
000)  , 
DISP  = 
000), 
OISP  = 
ono) , 
CL=SE 
ABEL  = 
OISP  = 
0  00)  , 


)  ** 
(  NEW 
SPAC 
(  NEW 
SPAC 
(NEW 
SPAC 
(NEW 
SPAC 
(MEW 
SPAC 
(MEW 
SPAC 
(NEW 
SPAC 
(NEW 
SPAC 
(MEW 
SPAC 
(  MEW 
SPAC 
(NEW 
SPAC 
R=DI 
(,  ,, 
(MEW 
SPAC 


,OEL 

E=(C 

,OEL 

E=(C 

,OEL 

E  =  (C 

,DEL 

E=(C 

,OEL 

E=(C 

,OEL 

E=(C 

,OEL 

E  =  (C 

,OEL 

E=(C 

,OEL 

E=(C 

,OEL 

E=(C 

,OEL 

E  =  (C 

SK01 

IN) 

,DEL 

E=(C 


-tt*  3^C  -nr*  3jC  nr  V*  -r^nr^ 

ET<= 
YL, 

ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 
ETE 
YL, 


ETE 
YL, 


5  , 
, 

5  , 
, 
5, 


5  , 


1)  ) 


**********    INSERT    ADIMA 
*********************** 

3LAN<    CARO    

BLAVK    CARO    

*************    DELIMITED 


****************************** 

DATA    HERE    *********** 
********************** 


(TW0    3LANK    CARDS    STOP    EXEC) 
CARO       (/*)     **************** 
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APPENDIX  D 
PSAP1   LISTING 
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